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The J1-J2 square lattice Heisenberg model with spin S = 1/2 has three phases with long-range 
magnetic order and two unconventionally ordered phases depending on the ratio of exchange con- 
stants. It describes a number of recently found layered vanadium oxide compounds. A simple 
means of investigating the ground state is the study of the magnetization curve and high-field sus- 
ceptibility. We discuss these quantities by using the spin- wave theory and the exact diagonalization 
in the whole J1-J2 plane. We compare both results and find good overall agreement in the sectors 
of the phase diagram with magnetic order. Close to the disordered regions the magnetization curve 
shows strong deviations from the classical linear behaviour caused by large quantum fluctuations 
and spin-wave approximation breaks down. On the FM side (Ji < 0) where one approaches the 
quantum gapless spin nematic ground state this region is surprisingly large. We find that inclusion 
of second order spin- wave corrections does not lead to fundamental improvement. Quantum cor- 
rections to the tilting angle of the ordered moments are also calculated. They may have both signs, 
contrary to the always negative first order quantum corrections to the magnetization. Finally we 
investigate the effect of the interlayer coupling and find that the quasi-2D picture remains valid up 
to \J±/Ji \ ~ 0.3. 



PACS numbers: 75.10.Jm, 75.30.Ds, 75.60.Ej 

I. INTRODUCTION 

The search for a quantum spin-liquid — an insulating 
magnet vifith a gapless ground state which breaks nei- 
ther lattice nor spin symmetries — has focused largely 
on spin-1/2 two-dimensional quantum antiferromagnets 
(2DQAF's). In practice, however, most two-dimensional 
AF spin-1/2 Heisenberg models exhibit either Neel order, 
or crystals of short-ranged singlet bonds and a finite gap 
to spin excitationsi. In a few cases, gapless hidden order 
states with nematic character arise^, but among "realis- 
tic" models perhaps only the S=l/2 Heisenberg model 
on a Kagome lattice remains a serious candidate for a 
spin-liquid description (see, e.gi^ and references therein). 

Perhaps the best studied example of a 2DQAF is the 
spin-1/2 J1-J2 Heisenberg model, which demonstrates 
a quantum phase transition from Neel to valance bond 
solid as a function of the control parameter J2/J1— . A 
number of layered vanadium compounds have recently 
been synthesized which are well described by this model. 
They are of the type Li2VOX04 {X = Si, Ge)^-^-^ and 
AA'VO(P04)2 (A, A' = Pb, Zn, Sr, Bei)2^ consisting of 
vandium oxide pyramid layers containing ions with 
spin 5=1/2. 

In Refs. lTollTll an extensive analysis of the J1-J2 model, 
also for finite magnetic field has been given in order to 
understand physical properties of the above compounds. 
Both the numerical exact diagonalization (ED) Lanczos 
method for finite clusters and the analytical spin-wave 
analysis have been employed. The behaviour of the sat- 
uration field as a function of the frustration angle has 
been studied as a further means of diagnosis of J1-J2 



compounds. It was found that close to the disordered 
regime with Ji < it is determined by an instability of 
two spin excitations, indicating indeed that in this regime 
the ground state may be of a spin-nematic typcj^i This 
leads us to expect that the magnetization itself should 
also be anomalous in this regime with a large effect from 
quantum fluctuations. So far quantum corrections to the 
magnetization curve in the spin-wave theory have only 
been considered for the nonfrustrated square-lattice an- 
tiferromagnet (J2=0)f^ 

In this work we give a systematic investigation of the 
magnetization and the high-fleld susceptibility for the 
general 2D square lattice J1-J2 model. Our goal is to 
investigate how the quantum corrections on these quan- 
tities depend on the degree of frustration, especially close 
to the disordered phases. The effect of interplane cou- 
pling will be also considered. 

The low scale of exchange interactions which are of 
order 10 K in the Li2VOX04 and and AA'YO(P04)2 
vanadates mean that magnetization measurements are 
relatively easy to perform. It is therefore hoped that the 
analysis with the theory developed here gives an addi- 
tional criterion to determine the frustration ratio J2/J1 
in a specific compound. In Sec.|II]we briefly introduce the 
model. Its high fleld properties like magnetization and 
susceptibility are investigated in Sec. IIIII They are ob- 
tained from ED Lanczos calculations and first as well as 
second order spin- wave theory and a comparison is given. 
We also evaluate the contributions of interlayer coupling 
in Sec. IIVI and in Sec.|V]we finally give a discussion and 
conclusion. 
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II. THE J1-J2 MODEL AND ITS PHASES 

The 2D square lattice spin-1/2 Heisenberg model in an 
external magnetic field H is given by 

n = JiY, s^-s, + J2 ^ s,.s, -hY,s: , (1) 

where Ji and J2 are the two exchange constants between 
the first and the second neighbors on a square lattice, 
respectively. As in Refs. [lOlilll the exchange parameters 
are defined per exchange bond. Furthermore we use the 
convention h — gfisH (g — gyromagnetic ratio, /i^ = 
Bohr magneton). The phases in zero field are best char- 
acterized by introducing equivalent parameters 

Jc = {J^ + Jl)K <t> = tan-i(J2/Ji) , (2) 

or j = tan (p = J2I Ji- The angle </> determines the extent 
of magnetic frustration in the model. 

Three classical magnetic ground states are possible de- 
pending on 0; namely ferromagnet (FM), Neel antiferro- 
magnet and coUinear antiferromagnet (CAF). They have 
been extensively discussed in Ref. T^. The effect of ex- 
change frustration leading to enhanced quantum fluctua- 
tions is strongest at the classical phase boundaries where 
the CAF phase joins the NAF or FM phase (see insets 
of Fig. [1]) . In fact in these regions they are believed to 
destroy long-range magnetic order and establish two new 
partially ordered states, namely a columnar dimer state 
with a spin gap at CAF/NAF boundaryi^iHiiiOSiie ^ 
gapless spin nematic state at the CAF/FM boundary-- . 
It can be seen already within spin wave approximation 
that the magnetic order breaks down in this regime since 
the sublattice moment reduction due to quantum fluctu- 
ations diverges close to the two boundary regionai^i. 



III. HIGH FIELD PROPERTIES OF THE J1-J2 
MODEL 

The determination of the frustration ratio Jij J\ or an- 
gle <p = tan^^(J2/ Ji) is of foremost importance to char- 
acterize a given square lattice magnetic compound such 
as the Si04 and PO4 vanadates mentioned in the intro- 
duction. The available experimental methods to reach 
this objective have been extensively discussed in Ref. 
[lOl . However, the thermodynamic zero-field methods like 
evaluating the heat capacity and the magnetic suscepti- 
bility are too ambiguous to locate in the NAF or CAF 
sector of the phase diagram. Additional information may 
be obtained from investigating high-field properties'^, for 
example from saturation fields as, e.g., determined from 
the magnetocaloric effect. 

The analysis of the high field magnetization itself as 
function of frustration angle is also promising. From 
the simple n.n. Heisenberg AF (J2 = or = 0) it is 
known from analytical work'^ that deviations from the 



classical linear magnetization curve due to quantum fiuc- 
tuations are to be expected. This is also concluded from 
numerical calculations A systematic study of magne- 
tization curves for the J\-Ji model is however lacking. 
Since the method is of experimental importance due to 
its relative simplicity we consider it worthwile to inves- 
tigate this problem in detail for the antiferromagnetic 
phases of the J1-J2 model. For this purpose we will use 
both analytical spin wave methods similar to Ref. [H and 
numerical Lanczos methods for exact diagonalization of 
finite clusters. Our notation will be close to the one used 
in Ref. [nl 



A. Magnetization from numerical {T = 0) Lanczos 
results for J1-J2 clusters 

We have numerically diagonalized the Hamiltonian of 
the J1-J2 model for cluster sizes of 16, 20, and 24 sites 
applying the finite-temperature Lanczos method as de- 
scribed in Ref. 11^ and references cited therein. In the 
limit T — > 0, our results are identical to the standard 
zero-temperature implementation of the algorithm, i.e., 
the evaluation of the partition function reduces to the 
determination of ground-state expectation values. The 
clusters of size 16 and 20 are (regular and tilted) squares, 
the 24-site cluster is a rectangle. All three clusters tile 
the infinite lattice such that, with periodic boundary con- 
ditions, compatibility with the three classically ordered 
ground states is preserved. 

The zero-temperature field dependence of the magneti- 
zation m = (S"^) has been calculated for the whole phase 
diagram except the FM region, where the saturation field 
hs vanishes. Following Ref. [20, Fig. [1] shows the normal- 
ized magnetization m/mgat (wsat = S) for selected values 
of the frustration angle (fi in the classically ordered anti- 
ferromagnetic phases and in the two disordered regimes 
of the phase diagram. The magnetic field is normalized 
to the saturation field determined by exact diagonaliza- 
tion. This is identical to the classical saturation field for 
positive (antiferromagnetic) Ji. For negative (ferromag- 
netic) Ji, a AS" — 2 two-magnon instability determines 
the saturation field for the finite-size systems considered 
here, which occurs at slightly higher field values than the 
one-magnon instabilityfii 

The field dependence of the magnetization at T = 
for a finite-size system is a sequence of finite steps. The 
solid line in the left panel of Fig. [T] shows the magnetiza- 
tion curve of the 24-site cluster for (p — —0.21 tt, which 
is in the Neel phase, the dashed line shows the same for 
(f) = 0.75 TT in the coUinear phase. In the right panel, 
the solid line shows the field dependence for ip/ir = 0.17 
(columnar dimer phase), the dashed line for (p/n = 0.84 
(spin-nematic phase). The symbols in the plots denote 
the midpoints of the horizontal and vertical line segments 
of the magnetization steps; circles label the 24-site data, 
diamonds the 20-site data, and squares the 16-site data. 
(For the smaller cluster sizes, the step functions are not 
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FIG. 1: The normalized magnetization obtained from the Bonner-Fisher construction as a function of an applied magnetic 
field H in units of the saturation field H^at for diff'erent frustration angles (for Ji > Hsat = g/J-Bhs). Data from the 24-site 
cluster are denoted by circles, 20-site data by diamonds, and 16-site data by squares. Additionally, the zero-temperature 
magnetization steps are plotted for the 24-site clusters; the solid line corresponds to the solid symbols, the dashed line to the 
open symbols. The insets in the plots show the positions of the frustration angles (p/n in the classical phase diagram. The 
frustration angle is counted positive above (0 < 4)/tt < 1) and negative below ( — 1 < ^/tt < 0) the x-axis in the inset diagram. 
For reference, the classical magnetization curve is plotted as a dotted line in both plots. Left: Classically ordered phases. Neel 
phase (NAF), (/> = — 0.21 tt (solid symbols) and colhnear phase (CAF), </> = 0.75 vr (open symbols). Right: Disordered regimes. 
Columnar-dimer phase (crossover from NAF to CAF), = 0.17 tt (solid symbols) and spin-nematic phase (crossover from CAF 
to FM), (f) = 0.84 TT (open symbols). 



shown.) Note that due to the AS' = 2 steps in the mag- 
netization for (j) > tt/2 (ferromagnetic Ji) there exist 
only half as much data points than for < 7r/2. Ac- 
cording to Bonner and Fisher (Ref. [20l). these midpoints 
should yield, at least for J2 = 0, a good approximation of 
the magnetization curve of the square-lattice Heisenberg 
model in the thermodynamic limit. 

The angles c^/tt = -0.21 and 0.75 for the NAF and 
CAF in the left plot are chosen such that they corre- 
spond to the experimental findings for the compound 
SrZnVO(P04)2.^° For these values, the midpoints of the 
magnetization steps form a smooth function. This can 
generally be observed for any frustration angle located 
inside the magnetically ordered regimes, be it NAF or 
CAF. In contrast, for (p/n = 0.17 (solid symbols and 
line in the right plot) and (p/n — 0.85 (open symbols 
and dashed line) , the data points scatter much more and 
do not give rise to a smooth field dependence. Further- 
more, in the columnar dimer phase ((/'/tt = 0.17) the 
well-stablished half-magnetization plateau that appears 
at m = imsaf^"^- The data in the right panel of Fig. [T] 
indicate possibility of yet another plateau at m = jnisnt, 
though a careful examination of fintie-size effects is nec- 
essary to make a final conclusion. 

In the whole phase diagram, quantum effects lead to 
negative corrections: the exact magnetization curve m{h) 
lies always below the corresponding classical value, which 
is a consequence of the lowering of the ground state en- 



ergy of a quantum antiferromagnet compared to its clas- 
sical counterpart 

B. Spin wave excitations in an external magnetic 
field 

A standard Holstein-Primakoff approximation of 
Eq. ([T]) and a subsequent Bogoliubov transformation 
leads to the harmonic spin wave Hamiltonianii 

n = NEo + NE,p + J2 eAk(/j)alk«Ak - (3) 

Ak 

where a^^. are magnon operators that obey bosonic com- 
mutation rules. The £\.k{h) denote the spin wave disper- 
sion of branch A = ± as defined in the appropriate NAF 
or CAF magnetic Brillouin zone (BZ). It is given by 

e±k(M = V[^kW±Ck(/i)]2-Sk(/i)2 . (4) 

Here ^k(^) is the intra- and B]^{h), C^ih) are the inter- 
sublattice couplings given below. The ground state en- 
ergy is composed of a classical part (Eq) obtained from 
the mean field approximation to Eq. ([T]) and a part due 
to zero point fluctuations of spins {Ezp). The former is 
given by 

Eo = -/i(S,|)-|-a||(S||)2-a^(Si)2 , (5) 
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FIG. 2: Magnetization curves ^i/ ^lb = gm (= m/S) for various frustration angles in the antiferromagnetic or disordered sectors 
with an offset of 0.2 applied. Symbols are obtained from Lanczos magnetization data for N = 16 (squares), 20 (diamonds), 
24 (dots, circles) size clusters using the Bonner- Fisher construction^^. Lines are obtained from first (full) or second (dotted) 
order spin wave calculations of Sects. UlI CI and [ III Fl Angles (p/n = 0.75, -0.21 correspond to the possible CAF or NAF values 
of the Sr compound. Magnetization curves strongly differ in the extent of nonlinear deviation from the classical curve which 
corresponds to = -0.5. Furthermore i^/tt = 0.41, -0.17 are values which are deeply within the CAF or NAF regions and 
overall agreement of spin wave and Lanczos calculations is good. It is less so on the CAF side where the magnetization changes 
in steps of ASz—2 leading to a larger finite size scattering. The values (p/n = 0.75, 0.17 correspond to regions close to or 
within the nonmagnetic sectors. Close to the CAF/FM boundary the first order spin wave results overemphasize the nonlinear 
behaviour and become unstable at very low fields. Close to the CAF/NAF boundary the numerical data exhibit a plateau at 
m/S — 1-l/ij.b = 0.5 which would require a separate analysis. The plateau was first reported in Ref . [2^ The inset shows the 
position of plotted (p values in the phase diagram. Second order spin wave results are discussed in Sect. IIII Fl 
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where {Si\) = Scos^, {S±) = S'sin% and 6*^/2 is 
the classical field induced canting angle of AF moments 
counted from the z-axis which is chosen parallel to the 
applied field h (here ||= z and _L= x,y). The coefficients 
in Eq. ([5]) are given by 



ail 
a± 



:(Ji-J2 



(NAF) 
(CAF) 



(6) 



where 2: = 4 is the coordination number. The classical 
canting angle 0c/2 of moments is obtained minimizing Eq 
which leads to 



h 

h. 



2S'(a|| + aj_) 



(7) 



Due to the symmetry properties ^k+q = ^k, ^k+q = 
— Bk and Ck+q = — Ck we have the identity e±(k+Q) — 
eip(k). Then, instead of summing over two spin wave 
branches in the magnetic BZ of NAF or CAF, we may 
restrict to one mode only, e.g. ek = e+k and sum over 
the whole paramagnetic BZ in the expression for E^p and 
similar ones. Using this convention the spin wave branch 
index A will be omitted in the following . Using Eqs. @ 
and ([9]) the spin wave energies may be written as 



ek(^) = 5'(ak + Ck)2(ak + CkCos6'c)2 , (12) 

for both NAF and CAF cases with ak and Ck given in 
Eq. (Uni). 



where hg is the classical saturation field. For h> hg the 
moments are fully polarized, i.e., all are ferromagnetically 
aligned parallel to h. Explicitly we have hg = 2zSJi 
(NAF) and hg = zS{Ji + 2 J2) (CAF). 

The zero point energy due to quantum fluctuations is 
obtained asii 



(8) 



Ak 



For the CAF and NAF phases it is always negative and 
vanishes in the FM or fully polarized (h = hg) phase 
where the ground state and spin wave states are exact 
eigenstates with a dispersion e^k = ^k- The sublattice 
couplings which determine the spin wave dispersion in 
the canted state were derived in Ref. [ll| and are given 
here for completeness. Defining 



C. Magnetization from first order spin wave 
quantum corrections 

The zero temperature magnetization is given by the 
field derivative of the total ground state energy: 



m 



mo + m,zp 
hg 2N^ 

k 



dh 
dh 



dh 



(13) 



where the first term is the linear classical part and the 
second one the (negative) correction due to quantum fluc- 
tuations included up to first order in 1/5. We can write 



S'flk 



Bk{h) = S'fokSin^ y , 
Cu{h) = 5ckCos2^ , 



(9) 



we have in the NAF and CAF phases corresponding to 
wave vectors Q — {tt,tt) and Q — (7r,0), respectively: 

flk = 4[Ji- J2(l-7k)] , (NAF) 
Ck = -fok = 4Ji7k , (10) 
flk = 2[2 J2 + Ji-/y] , (CAF) 

Ck = -fek = 2[Jl + 2J2jy]jx , 

where the geometric structure factors are defined by 



7k = ^(cosfc^ + cos ky) , 

7k = cos kx cos ky , 

7a; = COS kx, 7y = cos ky . 



(11) 



dh 



1 2Ck 



ek(/i) 
2S ( 
h-g'< 



[(Ak + Ck) - Bk] 



— COS — . (14) 

Ok + Ck COS Qc' 2 



This finally leads to a total magnetization, including the 
first order quantum corrections: 



TO 



1 1 

'hg'N 



E 



Ck 



ak + Ck 
' flk + Ck cos ( 



, (15) 



where on the r.h.s. the classical value of 6'c given by 
cos(0c/2) = hjhg has to be used. Because hg ~ 2S'(a|| + 
a_L) the second term in Eq. is formally a 1/5* correc- 
tion to the linear classical term toq = S{h/hg). 

Explicity, using Eqs. (jlOllip we have for the NAF case 
a magnetization depending on field strength and frustra- 
tion angle according to 
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FIG. 3: Left panel: Susceptibility as function of field for various ^/tt values in the NAF sector values correspond to curve 
sequence for h = 0.) Note the divergence close to ha for (\> approaching the sector of the stacked dimer phase. In this regime the 
zero field susceptibility tends to zero indicating the breakdown of the (1/S) approximation. Right panel: Zero-field normalized 
susceptibility Xn in the </!> > part of the NAF sector. At (/> = (see also left panel) the first order Xn (dashed line) obtained 
from Eqs. (I20|l and (|22|l is already reduced to half the classical value (equal to S). It decreases further with increasing and 
then becomes negative, indicating the instability of the NAF state. Inclusion of second order contributions in (1/S) stabilizes 
the positive value (full line, symbols) but they diverge on approaching the CAF/NAF boundary around (/j/tt ~ 0.15. 



NAF 



2SN 



1 +7k - j(l -7k) 
I- 7k cos 6*^ - j(l -7k) 



(16) 



where we used j — tan (j) = J2/ Ji- In this expression 
the {1/S) character of the quantum correction becomes 
manifest. For the simple NAF with j = we reproduce 
the result first given in Ref . il^ 



A similar but more complicated expression may be 
given for the CAF phase. Defining (5k = (1/2) (cos fc^, — 
cos ky) in addition to Eq. pT|) we obtain 



J 



CAF 



m ~ b — 
ft.. 



[o(7k + 5k) 



j(l + 7k) +7k 



j(l + 7k COS 9c) + 7k cos^ ^ 



5kSin2% 



(17) 



r 



Note that the special CAF case with Ji = {j = 00) 
is equivalent to the simple NAF case J2 ~ {j = 0) in 
Eq. p^ . This may be seen by applying a k-coordinate 
rotation by 7r/4 in Eq. p7|) . 

Naturally, the above expressions should be primarily 
valid deep inside the NAF and CAF regions where the 
staggered moments are large. Close to the boundaries 
quantum fluctuations grow and destroy the magnetic or- 
der, then corrections to m{h) starting from the ordered 
state and expanded in orders of 1/S are no longer appro- 
priate. 

The combined analytical spin wave and numerical 
Lanczos results for the magnetization are shown in Fig. [2l 
One may indeed see that the agreement of both is good 
deep inside the NAF (^/tt = -0.21,-0.17) and CAF 
{(j)/Tr = 0.41) sectors. On the other hand close to the 



classical phase boundaries {(f>/n — 0.17,0.75) discrep- 
ancies appear. For the former case our and previousS^ 
Lanczos results indicate appearance of the one-half mag- 
netization plateau with m = ^mgat close to h/h^ = 0.5. 
(Such a plateau has been also found in the equivalent 
classical model at finite temperatureSfSSi For the latter 
the magnetization curve becomes very nonlinear due to 
quantum fluctuations and m{h) becomes negative at low 
fields (upper dashed line) indicating the breakdown of 
the 1/S expansion. This means that the CAF state be- 
comes unstable and a new (spin nematic) order parame- 
ter will be realized in a finite sector around ^/tt — 0.85 
{J2/J1 ^ —0.5). Generally, one may say that quantum 
corrections leading to nonlinear magnetization will be 
considerable larger on the ferromagnetic ( Ji < 0) side of 
the CAF sector. This will be further discussed in Sec. IVl 



7 



D. First order quantum corrections to the 
magnetic susceptibility 

In the spin-wave approximation the magnetization 
curve for the simple AF exhibits a logarithmic singular- 
ity of the slope close to the saturation field hg as shown 
in Refs. 'l8"2l'. This should be more easily visible in the 
high field susceptibility. Furthermore for cj) approach- 
ing a classical phase boundary the low field susceptibil- 
ity has to vanish, which signifies the instability of the 
order parameter. For these reasons we found it useful to 
study the magnetic susceptibihty x{h, (j>) = dm{h, (j))/dh 
as function of the frustration angle in addition to the 
magnetization. From Eq. (|13p wc obtain 



X{h, (t>) =Xo + Xzp 



S_ J_v- 



(18) 



From Eq. (|12p and using {dA\s_/dh) = we arrive at 



\ dh ) 



1 / de^{h) 
e^{h)\ dh 



Inserting this in Eq. (fT8|) and using Eq. (fT4|) leads to 
a general expression of the normalized susceptibility 
Xn{h) — hsx(h) according to 



Xn{h,4>) = S 



l--(Ax(")(/i,0)-Axi'')(/i,<^)) 



(20) 



where the terms ~ l/S* in brackets are the quantum cor- 
rections of to the constant classical value x° = S. 
We obtain from Eqs. (|18ll9p : 



^X^\h,c^) 
Axl^HKcj^) 



Qk + Ck 

Ctk + Ck cos ( 



(21) 



hAhJ N 



(flk + Ck)^ 

(ok + Ckcos 9c)i 



(19) For the NAF case this may be explicitly evaluated as 



^xi^Hh, 



1 

2N 



-Y-y 



1 + 7k ~ j(l -7k) 
f 7k cos 6*^ - j{l - 7k)' 

,2 (l+7k-i(l-7k))^ 
(1 7k cos 6'c - j(l -7k))^ 



(22) 



Obviously only the first part A^n contributes to the 
(l/S) corrections of the zero- field susceptibility. It may 
also be directly obtained by differentiation of m(h) in 
Eq. PB)) . For the CAF case similar expressions for 

AxrT''''' {h, (f)) may be derived by making analogous sub- 
stitutions in the integrals and their prefactors as done in 

Eqs. (HH) and ini). 

The typical field dependence of the susceptibility is 
shown in Fig. H For cp/n = -0.49 close to the NAF/FM 
boundary one obtains nearly the classical constant value 
Xn — S because in the FM sector quantum fluctuations 
are not present, they are gradually turned on when Ji 
becomes positive and (j) moves into the NAF sector. This 
can be clearly seen from the various curves in Fig. [3] (left). 
For an angle ^/tt = 0.1 frustration becomes large and 
quantum fluctuations are close to destroying the NAF 
order. Accordingly the zero-field susceptibility is close 
to becoming negative where the spin- wave theory breaks 
down. One also notes the upturn in the susceptibility just 
below the critical field coming from the logarithmic singu- 
larity of the magnetization. 1*'--'^ The singularity becomes 
stronger when the strongly frustrated point j — 1/2 is 
approachedi^i 



E. Quantum corrections of the moment canting 
angle 

The angle 9c between canted AF ordered moments has 
sofar been given in classical approximation from minimiz- 
ing only Eq(0c, h). Its quantum corrections of first order 
in (l/S) may be computed by minimizing the total en- 
ergy Eo{9c, h) + Ezp{9c, h), keeping 9c also as variable in 
the zero point energy. For the unfrustrated Neel AF this 
correction is small, however, as we shall see in the follow- 
ing it may be of considerable size close to the classical 
boundaries of NAF and CAF phases where frustration 
effects are large. 

The equilibrium condition including quantum correc- 
tions is given by 

d0^+2N^^d9c d9c'~ ^ ' 

k 

Now 9c is to be treated as a variable present in ek and ^k- 
The form of in Eqs. (|9I10[) has already the classical 
angle of Eq. ([7]) substituted. Its general form is Ak = Sa^ 
with 

Q 

flk = -4J2(1 - 7k) - 4Ji cos 9c + h cos ^ , (24) 



8 




0.2 0.4 0.6 0, 
h/h 



1 0.2 0.4 0.6 0.8 1 

h/h 



FIG. 4: Comparison of (1/S) quantum corrected canting angle 6'^ (left) and magnetization (right) for two values in the NAF 
sector. It is seen that first order quantum corrections always reduce the magnetization with respect to the classical value 
according to Eq. (|15[) . For cos(^^/2) however the corrections may have both signs. They are small and negative for (j) deep in 
the NAF sector while they are positive closer to the NAF/CAF boundary. 



for the NAF case. Likewise for CAF sector one obtains Evaluating the derivatives of Ck and with respect to 

Oc and solving the equilibrium Eq. (j23p we obtain the 
'Ik = ^2Ji(l — ^fy) canting angle 9'^ renormalized by quantum fluctuations: 



-2( Ji + 2 J2) cos 9c + h cos ■ 



(25) 



--y\ 

hsN ^\ 



,(1) 



(flk + Ck cos^ + Ck(ak + Ck) 
(flk + Ck) 2 (ak + Ck cos 9c) 2 



,(1) 



(26) 



r 



Here cos — h/hg is the classical canting angle and 

a\^^ = —{hs/2S) is the coefficient of the second term in 
Eqs. (|24l25p . Again, because of hg = 2S{a\\ + a±) the 
term ^ l/hg in Eq. ((26)) has to be considered as a (1/S) 
correction to the canting angle. Therefore in the smn on 



the r.h.s. the classical values for 9c has to be used which 
leads again to the intra-sublattice coupling Ok as given 
by Eq. (Uni). 

The general solution in Eq. ([26]) is valid for both AF 
phases. In the NAF case it leads to the explicit result 



J 



cos — = cos 



2SN ^ 

k 



7^ 



- 7k sin^ |l - 1 + j(l - 7k)(l - 7k) 



(1 +7k - j(l - 7k)) 2 (1 + 7k cos 61c -j(l -7k))' 



(27) 



r 



For the simple AF (j = 0) one recovers the expres- 
sion given in Ref. [isl . A similar explicit expression for 
the CAF case may be derived but it is too unwieldy to 
be given here 

Eq. ([26]) may be used as well 

The field dependence of renormalized moment cant 



For the numerical calculation of cos 



ing given by cos in comparison with the normalized 
magnetization m(h)/S is shown in Fig. H) Generally the 
quantum corrections to the canting angle are quite small 
deep inside the AF sectors. As for the magnetization they 
become larger when approaching a phase boundary. In- 
terestingly however, they may have different signs for the 
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former (left panel) while they must always be negative for 

the latter (right panel) . The positive correction to cos ^ 
appears in the vicinity of the NAF/CAF boundary. 

The quantum corrections to magnetization and cant- 
ing angle may be used to obtain the correction to the mo- 
ment size. Classically we have S = mo/ cos ^ . If we use 
a similar definition including the quantum corrections we 

have (S) = m/ cos ^ for the renormalized moment. The 
change of moment size defined hy SS = S — (S) is then 

given hy SS = S — [raj cos ^). Note however that in this 
relation the quantum correction to S formally contains 

effects of arbitrary order in (1/S) even if m and cos 
are only corrected in order (1/S). Close to the nonmag- 
netic regions when m aproaches zero SS/S becomes unity, 
i.e., the staggered moment is destroyed by the quantum 
fluctuations. 



F. Second-order spin-wave results for the 
magnetization and susceptibility 

The linear spin wave theory used in the previous sec- 
tions includes the leading contribution to the ground 
state energy 0{S), see Eq. ([3]). The next order contri- 
butions to -Eg.s. come (i) from the canting angle renornial- 
ization, (ii) from a Hartree-Fock decoupling of the quar- 
tic terms, and (iii) from the cubic terms. The detailed 
derivation in the case of nonfrustrated square lattice AF 
is given in Ref. Ts!. Here, we present only the final result 
for the Neel phase of the Ji- J2 model. The ground state 
energy per site in second order spin-wave theory consists 
of the classical part Eq, first order corrections (zero point 
fluctuations) Ezp and second order in 1/S' corrections. In 
total it may be written as 



I6J1 2N 



9.N ^ 



2 Ji [(n - 8xf + ni(ni - S)\ + 2 J2 {{n - n^f + 82{S2 - S)\ 



+ 4Ji cos2(6lc/2) [n ~ + nx) ~ 25ini 



Ji sin^ Qc 



(28) 



F(k,q)2 



q + £k+q-Q 



r 



Ezp while the 



Here the first line is equal to the sum of Eq 
second and third line constitute the second order in 1/5 
spin-wave contribution. In this term the six constants 
are given by two-dimensional momentum integrals 



-E 

N ^ 



Ak + Ck - ek 



N ^ 2£k 

k 

"1(2) = ^E'^k(7k) 
k 



^k + Ck 

2£k 



(29) 



^ = ^E|^' '^i(2) = ^E^'^(7k)|^, 

k k 

while the expression for the cubic vertex i^(k, q) is given 
by Eq. (25) in Ref. [H. The magnetization curve is ob- 
tained by numerical differentiation: vn = —dEg ^X^)/'^^- 
The effect of the second order contributions in the mag- 
netization may be seen in Fig. [5] and also in Fig. [5] in the 
comparison to the ED Lanczos results. In general the sec- 
ond order spin- wave corrections to the first-order in {1/S) 
results are positve for small fields and negative for large 
fields (Fig. [5]). As long as (j> is within the stable antifer- 
romagnetic sectors the second order corrections are quite 
small. For example when (p/ir — -0.21, -0.17, 0.41 first- 
and second-order results are almost indistinguishable to 
the eye and agree with Lanczos results (Fig. [J). When 
4> approaches the strongly frustrated phase boundaries 
first-order approximation breaks down as witnessed by 



the magnetization becoming negative. Figs. [2] and [S] The 
second-order corrections remedy this situation. However, 
close to the boundaries they become very large leading 
to a very anomalous low field second order magnetiza- 
tion which deviates strongly from the Lanczos results 
(Fig. [2]) . Since the deviations between first and second 
order curves are strongly enhanced close to the classical 
phase boundaries it is clear that spin wave expansion no 
longer converges. In Fig. [2] this is obvious for (/)/7r = 
0.17 at the NAF/CAF boundary and even more for (p/n 
— 0.75, already well in advance of the CAF/FM bound- 
ary at (j>/Tr — 0.85. In fact the first order result (long- 
dashed) agrees better with the Lanczos results (full sym- 
bols) plotted for comparison than does the second order 
curve (short-dashed). 



Finally, we mention that from Eq. ([29)1 the second 
order in (1/5) corrections to the zero-field suscepti- 
bility may be calculated in addition to the first order 
expression given in Eqs. (|20|) and (|22p. The comparison 
of first and second order results for the NAF case is 
given in Fig. [3] (right panel). Although the second 
order contributions repair the negative instability of the 
first order susceptibility at (j> — O.IItt the second order 
result itself diverges when one moves even closer to the 
classical phase boundary at ~ 0.157r. 
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FIG. 5: Compiled comparison of first (full lines) and second (dotted lines) order in (1/S) spin wave results for the magnetization 
for various frustration angles indicated in the inset. An offset of 0.2 has been applied. Curves in increasing order correspond 
to (p- values in counter-clockwise direction. Well inside the antiferromagnetic regimes {(j>/iT = 0.093, 0.25, 0.65) first and second 
order results show little difference. Close to the classical phase boundaries (^/tt — 0.135,0.16,0.75) first order results lead to 
instable (negative) low-field magnetization. The instability region below 4)/tt = 0.85 (CAF/FM) is much larger than around 
(j)/TT = 0.15 (NAF/CAF). Second order results are are always positive but show very anomalous low field magnetization contrary 
to the Lanczos results in Fig. [2] Note that the second order corrections in (1/S) to the first order curves are positive for low 
fields and negative for large fields. 



IV. EFFECT OF INTERLAYER EXCHANGE 
COUPLING 



In real magnetic systems other interactions may play 
a certain role besides the in-plane Heisenberg exchange. 
These include various anisotropics as well as a three- 
dimensional coupling. Below we shall consider modifi- 
cations of the above formulas produced by interlayer ex- 



change coupling: 



H± — Jj 



E 



(30) 



for nearest-neighbor spins in the direction perpendicular 
to the layers. (Note the different convention for field di- 
rection in Sect. IIII Bl) . Here we assume the simple stack- 
ing of the layers. The obtained results can be easily ex- 
tended to other cases as well. 
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For ferromagnetic exchange Jx < 0, the saturation order spin- wave result for the magnetization in the NAF 

field for NAF is still given by hg = 8SJi, while for an- case which includes the interplanar coupling J± is given 

tiferromagnetic exchange > it is given by hg = by 
8S'Ji(l + ^j±). Here, we define j± = J±/Ji. The first- 



Jj_ < : m 
Jj_ > : m 



= S 



4 



1 - 



1 - 



2SN ^ 

k 



7k 



(r 



1 + 7k - j(l - 7k) + 5 1 i_L I (1 - cos fcz 



1/21 



2SN 



+ 7k cos 
V- ??k / 1_ 

4-i + i7lVi + i 



j(l-7k) + 5l,7±|(l-cosfc,) 
jXl-7k) \i/2 



(31) 



1 + + '7k 



+ t J_L ^ 1 + + ?7k cos ( 



j(l -7k) 



Here ?7k = 7k + \ cos is the 3D structure factor and 
summation is now extended over a 3D Brillouin zone. 

In the case of CAP order the first-order spin-wave the- 
ory for a ferromagnetic interlayer exchange Jj_ < leads 



to hg = 45Ji(l + 2j) and for antiferromagnetic case J± > 
the saturation field is given by hg = 4:SJi{l + 2j + 4:j±). 
The first order spin wave result for the magnetization in 
the 3D case is obtained as 



Jj_ < : m 
J_L < : m 







4 


1 


hs 








4 


1 


hg 





1 - 



1 



2SN 

1 

2SN 



E 

k 

E 



i7k + ^7x 
i + l 



(id 



i(i + 7k) + 7k + ||i±|(i 



1/2- 



i(l +7kCos6'c) + 5(7j/ + cos6ic7x) + \\3i.\{^ - 
j(l + 7k)+7k + 5i±(l- 



5i-L7^ 



(.7(1 



7.)^ 

H7.) 



(32) 



1/21 



3 + 5(1 + J-l) Vj(l +7kCos6'c) + 5(7y + cos6ic7x) + 5i±(l + 7z cos 6*0) > 



Here we used the convention 7^ — cosfc^ (« = x, y, z). 

We have verified that for well ordered NAF and CAF 
phases the effect of interlayer coupling is hardly visible 
up to ^ 0.3Ji^2- This means that the application 
of our first order spin wave results to real compounds 
does not critically depend on a small interlayer coupling 
J_L as long as the quantum antiferromagnet is in a well- 
ordered phase. Closer to the phase boundaries Jx has 
the the effect of stabilizing the ordered phases, especially 
for J_L < 0. 



V. DISCUSSION AND CONCLUSION 

In this work we have explored the quantum corrections 
to the magnetization, the uniform susceptibility and the 
canting angle in the first- and the second-order spin-wave 
approximation and by using exact diagonalization of fi- 
nite clusters. The deviations from classical behaviour 
were foimd to be pronounced close to the strongly frus- 
trated regions where the classical phases meet. Indeed 
we have shown that linear spin- wave theory breaks down 
in these regions and second-order corrections do not fun- 
damentally change this observation. The latter are small 
well within the AF regions and have positive sign for 
small and negative sign for large fields. Although they 



prevent the instability of the linear spin- wave theory for 
small fields they become very anomalous close to the 
boundaries and one cannot expect that the 1/S expan- 
sion converges in the nonmagnetic region. This is intu- 
itively clear since the spin wave expansion in the region 
of the classical boundaries starts from the wrong, i.e., 
magnetically ordered ground state. 

A striking feature of these results is that the deviations 
from classical behaviour in the magnetization curve m{h) 
and the breakdown of the spin-wave expansion are most 
pronounced on the ferromagnetic side (Ji < 0), as il- 
lustrated in Figs. 2 and 5. One can understand this by 
considering the asymptotic form of the first-order cor- 
rections to the staggered moments 5S and the uniform 
susceptibility 5x in the CAF phase as j — > ±1/2: 



5S 



5x 



E 



1 



I 

J + ^kl + \5j\kl 



(33) 



Here, 5j = j ± \ is the deviation from one of the two 
strongly frustrated points. For the CAF/NAF boundary 

the diverging corrections arc the same for the sublattice 
magnetization and the susceptibililty: 5S,5x ~ ln|(5j|. 



12 




10 



^ 4 



1 1 1 1 

. FM 


1 1 1 1 1 1 

NAF 


1 1 1 1 1 1 

CAF 




1 

-M. 


- 






c 




- 

















o 
o 




1 1 1 1 


1 1 1 1 1 1 


a 

h 

1 1 1 1 1 1 


D 

1 





-1.0 -0.5 0.0 0.5 1.0 



FIG. 6: Left; Variation of the moment fi/ I-Lb ~ gm as function of (p for h/hs=0.58 close to half-saturation. Circles are obtained 
from the Bonner-Fisher plots based on numerical (T=0) Lanczos results. Full lines include the (1/S) quantum correction from 
spin wave excitations according to Eq. (|15|l with respect to the constant classical value ^J./ ^J.B ~ m/S — 0.58. At the FM/NAF 
boundary where the quantum fluctuations are small the magnetization starts at this value and then is continuously reduced as 
one approaches the strongly frustrated sector around NAF/CAF boundary. There the magnetization jumps to a slightly larger 
value and then breaks down close to the spin nematic region around CAF/FM boundary. Right: Susceptibility x(0)/x(0 = 0) 
for h = hj {cj}) normalized to the simple NAF value. Circles are ED Lanczos results (ksT/Jc=0.2) and full lines are T=0 spin 
wave results according to Eq. (f20)) . 



For the CAF /FM frustration point the susceptibihhy cor- 
rection acquires an additional diverging prefactor: dx ~ 
(In \6j\)/{6j), which indicates that a long-range magnetic 
order is destabilized in a much wider window of J2/ Ji for 
this sign of Ji. 

This is in accordance with the ED Lanczos results 
where a tendency to bound state formation of spin waves 
as indicated by the ASz — 2 steps in the magnetisaton of 
finite clusters is observed^ ^. This is evidence that around 
~ 0.857r or J2/J1 — —0.5 the ground state will be of 
the spin nematic type as proposed in Ref. ITtI . It may 
be viewed as a quantum gapless phase with a Goldstone 
mode describing the collective long range excitations of 
a nonlocal quadrupolar order parameter. A second or- 
der transition between and the CAF phase and this spin 
nematic is permitted by the symmetry of the order pa- 
rameter, and the smooth evolution of ED spectra across 
the transition suggest that a second order does in fact oc- 
cur. The pronounced quantum fluctuations seen in spin 
wave theory lend further support to this idea. 

This can be most clearly seen in Fig. [5] (left panel) 
where the magnetization is plotted as function of the 
frustration angle 0. The decrease in magnetization from 
the classical value m/ S = h/hc at around half saturation 
{h/hc = 0.58) characterizes the strength of quantum fluc- 
tuations. Their effect increases from zero at FM/NAF 
boundary to a maximum at NAF/CAF boundary where 
a discontinuous jump in the magnetization occurs. For 
(j)/Tr > 0.5 (the CAF regime with J2 < 0) the reduction of 
m/S due to quantum fluctuations rapidly becomes large 
and the expansion in 1/5 breaks down. This is also seen 



in Fig. El 

The exact numerical results (circles) obtained from the 
Bonner-Fisher plots (hence the steps) as in Fig. [T] also 
show the strong reduction of the magnetization close to 
classical CAF/FM boundary . This is a signature of the 
true spin-nematic quantum ground state which does not 
break time reversal symmetry. Therefore it has no first 
order (linear) coupling to the magnetic field, resulting in 
a small magnetization. The right panel of Fig. [6] shows 
the comparison of high-field susceptibihty a.t h = h~{(l)) 
from first order spin wave theory (full line) at T=0 with 
the ED Lanczos results at small but finite T, plotted as 
function ot the frustration angle 4>. The susceptibility in 
both cases is normalized to the pure NAF case (0=0). 
The overall agreement of cj) dependence is quite good be- 
cause due to high fields the effect of quantum fluctuations 
is suppressed. Again the deviations are strongest close 
to the classical CAF/FM boundary where the saturation 
field approaches zero and the spin wave approximation 
breaks down. 

The deviations from classical results and the break- 
down of spin wave approximation is much less severe on 
the CAF/NAF phase boundary at the AF side (Ji > 0). 
The magnetization behaves significantly less singular in 
this region. As mentioned before a discontinuous jump in 
m(0) at ~ 0.157r occurs around half saturation (Fig.lH]). 
Then the quantum phase transition to the presumably 
stacked spin dimer ground state may be expected to be 
a first order transition. Therefore the magnetic and non- 
magnetic phases on both sides of the boundary will cor- 
respond to stable local minima of the free energy and 
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fluctuations will not be very pronounced leading to a less 
singular magnetization behaviour. The first-order nature 
of transition between CAF and a nonmagnetic columnar 
dimer state was noticed in early numerical workfi^ while 
for the boundary with the NAF, a first-order scenario 
was put forward only recently.— 

Finally, we comment on experimental data for the 
high field magnetization of J1-J2 compounds. There 
are no published data for those systems mentioned in 
the introduction. However recently a new compound 
(CuBr)LaNb207, with a perovskite/metal halide inter- 
growth structure was synthesized^ . This is a spin-1/2 
magnetic insulator, and is reported to exhibit quasi-2D 
magnetic behavior arising from the CuBr - square lat- 
tice planes. It shows CAF order at a relatively large 
Tn = 32 K. 

It has been suggested that the magnetism of 
(CuBr)LaNb207 can be described by the square lat- 
tice J1-J2 model with a frustration angle <j) — 0.737r 
(J2/J1 = —1-1), which puts this compound closer to the 
strongly frustrated region of the spin nematic phase than 
any other compound reported sofar. Our spin wave calcu- 
lations predict a pronounced nonlinear magnetization for 
this parameter set, particularly at low fields. However, 



the experimental results in Ref. '25' show only a mod- 
est curvature, at high fields. We therefore conclude that, 
while these materials likely do possess competing FM and 
AF interactions, they probably cannot be described by a 
simple square lattice J1-J2 model. 

In conclusion, it is abundantly clear that classically 
disordered phases of the square lattice model need 

and deserve an analysis which goes beyond the spin wave 
theory presented here. For further progress to be made 
in understanding the "hidden order" phases an approxi- 
mative treatment of the broken order parameter and its 
low lying excitations is necessary. Since both the dimer 
and nematic phases are bond centered, the bond-operator 
method might be a useful choice. 
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